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Abstract 

We analyze a class of supersymmetric models first introduced by Arkani-Hamed et al 
and Borzumati et al in which the light neutrino masses result from higher- dimensional 
supersymmetry-breaking terms in the MSSM super- and Kahler-potentials. The mech- 
anism is closely related to the Giudice-Masiero mechanism for the MSSM \x parameter, 
and leads to TeV-scale right-handed neutrino and sneutrino states, that are in principle 
accessible to direct experimental study. The dominant contribution to the light neutrino 
(Majorana) mass matrix is a one-loop term induced by a lepton-number violating S-term 
for the sneutrino states that is naturally present. We focus upon the simplification and 
analysis of the flavour structure of this general class of models, finding that simple and 
novel origins for the light neutrino mass matrix are possible. We find that a subdomi- 
nant tree-level 'see-saw' contribution may lead to interesting perturbations of the leading 
one-loop-induced flavour structure, possibly generating the small ratio Amg 0lar /Amj; tm 
dynamically. 



1 Introduction 



The case for the existence of small neutrino masses and associated physical neutrino 
mixing angles has enormously strengthened in recent years as a consequence of the now 
numerous experimental studies of atmospheric and solar neutrinos, and neutrinos from 
terrestrial sources. In particular, a recent global analysis [1] which incorporates data 
from CHOOZ [2], SNO [3], KamLAND [4] and Super-Kamiokande [5], leads to values 
for the mass difference squares, Am 2 , and three real mixing angles 8ij of the neutrinos 
in the 3a ranges, 1.4 x 10" 3 eV 2 < |Am 2 3 | < 3.7 x l(T 3 eV 2 and 0.36 < sin 2 # 23 < 0.67 
from the atmospheric data, and 5.4 x 10" 5 eV 2 < Am^ < 9.5 x 10~ 5 eV 2 and 0.23 < 
sin 2 02i < 0.39 from the solar data, while the remaining real mixing angle is bounded by 
sin 2 # 13 < 0.066. The phases that can appear in the neutrino mixing (MNS) matrix are 
currently unconstrained. 

The traditional and much studied explanation for these small masses is the see-saw 
mechanism[6]. This hypothesises the existence of two or more standard-model-singlet 
right-handed (rhd) neutrino states N iy with very large lepton- number violating Majorana 
masses, which couple to the weak-SU(2) lepton doublets Lj via a conventional Yukawa 
coupling XijL{NjH involving the electroweak Higgs. The Yukawa couplings are typically 
taken to be of size comparable to that of either the charged lepton Yukawas or the quark 
Yukawas, depending on the precise model. After the Higgs gains its vacuum expectation 
value, v, this Yukawa coupling leads to the Dirac mass matrix = \v. Finally, the 
light neutrino mass matrix obtained by integrating out the heavy rhd neutrinos is given 
by 

m u = -{m^) T M R l m^ (1) 

where Mr is the rhd Majorana mass matrix. If we take the heaviest Dirac mass to 
be of order either the r lepton or bottom quark mass, then agreement with the light 
neutrino mass inferred from the atmospheric neutrino anomaly requires roughly Mr ~ 
10 11 GeV, while if the heaviest Dirac mass is taken to be of order the top quark mass, 
as is commonly the case in explicit models, then M R ~ 10 15 GeV. The comparative 
proximity of this second scale to the inferred supersymmetric grand unified (GUT) scale 
Mgut — 2 x 10 16 GeV is interpreted as evidence in favour of the see-saw mechanism, as 
is the fact that rhd SU(3)xSU(2)xU(l) singlet neutrinos are naturally included in many 
GUT theories. 

Certainly therefore the see-saw mechanism is an attractive explanation of why the 
light neutrino masses are so small. However, it is not without its faults. In particular there 
is a tension between the strongly hierarchical nature of the observed Yukawa couplings in 
the quark and charged lepton sectors, and the essentially hierarchy-free masses implied 
by the Am 2 's. Moreover, both the #12 and #23 mixing angles are large while the angle #13 
is small which is in sharp contrast with the corresponding mixings in the quark sector 
which are all small. These problems can be solved in specific models, for example the Am 2 
values can be fitted by taking the spectrum of rhd neutrino masses to be hierarchical in 
such a way as to almost compensate for the hierarchical neutrino Yukawa couplings. But 
this has the price of introducing a wide range of rhd neutrino masses Mr ~ 10 10 — 10 15 
which then require explanation. 
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In addition, there is the worrying question of the testability of the see-saw mechanism. 
Given the large mass scale of the rhd neutrino states there is very little prospect of their 
ever being a direct test of the correctness of the see-saw mechanism. We are forced to fall 
back on indirect, and sadly not definitive tests. For example, the discovery of neutrino-less 
double beta decay would point to a Majorana nature for the light neutrinos as required 
by the see-saw mechanism, but as we shall see below, the see-saw mechanism is by no 
means the only way that such Majorana masses can be generated. 

In this letter we study a notably attractive alternative, advocated by Arkani-Hamed 
et al [7, 8], and Borzumati et al [9, 10] (for related work see [11]), that links the light neu- 
trino masses to TeV-scale supersymmetry breaking physics, and which has the significant 
virtue, compared to the see-saw mechanism, of being directly testable, at least in part, 
at the LHC and other proposed high-energy colliders. In overall structure our models are 
similar to those previously studied in Refs.[7, 8, 9, 10], in that the dominant contribution 
to the light neutrino masses arises from a one-loop diagram involving a supersymmetry 
breaking and lepton- number violating £>-term for the right handed sneutrinos, but, by an 
alteration of the model, we have been able to significantly simplify the way in which the 
flavour structure of the light neutrino mass matrix arises, and are thus able for the first 
time to study in detail some of the consequences of this very attractive class of models. 

2 Outline of the model 

Supersymmetric extensions of the standard model, such as the minimal supersymmet- 
ric standard model (MSSM), explain the origin of the weak scale in terms of the scale 
of the coefficients of the supersymmetry breaking soft operators that must be included 
within the MSSM. The usual assumption is that these susy-breaking operators arise 
from the (super)gravitational mediation of susy breaking that occurs primordially in 
some hidden sector at an intermediate scale rrij ~ 10 10 — 10 11 GeV, giving rise to soft 
susy breaking in the MSSM of order TeV ~ rrij/M. Here M is the reduced Planck 
mass M = M p \/\/fyK = 2 x 10 18 GeV. Infamously, there is one mass parameter in the 
MSSM, the jj,H u H d superpotential interaction, that naively appears to be independent of 
supersymmetry breaking. If this were indeed true then we would lose completely our un- 
derstanding of the origin of the weak scale in susy theories [12]. The realization of Giudice 
and Masiero [13] was that the potentially Planck-scale fj,H u Hd term in the superpoten- 
tial can be forbidden by a global symmetry, while an effective /z-term is generated from 
1/M-suppressed terms in the Kahler potential via supersymmetry breaking effects, thus 
naturally implying fi ~ TeV. 

As emphasized by Arkani-Hamed et al and Borzumati et al, the lesson of the Giudice- 
Masiero mechanism for neutrino masses in the context of susy theories is that SM-singlet 
operators, such as the rhd neutrino mass MrNN, or the neutrino Yukawa coupling 
XLNH U , might only appear to be renormalizable superpotential terms, but in fact may 
arise from 1/M-suppressed terms involving the supersymmetry breaking scale mi. 

Specifically, consider the usual MSSM Lagrangian to be supplemented by a set of 
superpotential and Kahler terms involving the rhd neutrino superfields iVj (i — 1, 2, 3 is 
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a generation index), and two Standard-Model-singlet chiral superfields X and Y which 
arise from the hidden sector. In general, the fields which communicate supersymmetry 
breaking to the neutrinos can either be flavour singlets or flavour non-singlets. Let us 
generically call the flavor non-singlet fields, X, (which therefore carry generation indices 
and the flavour singlet fields, Y, and for simplicity suppose that there is just one X 
field and one Y field. Then the model we wish to study is defined by the A^-dependent 
terms in the superpotential 

C™ = l d 2 6 [g^UN 3 E u + g'^L t N t H u + ...), (2) 

while the set of terms involving the rhd Ni fields in the Kahler potential are 

r ( y\ _ y\y . Y^YX- \ 

C% = J dH IhjjNiNi + h^NlNi + fo-^JV^- + . . . j . (3) 

The ellipses in Eqs.(2) and (3) stand either for terms involving the replacement of Y 
fields by Xij fields (with obvious changes to Ni flavour indices), or for terms higher order 
in the 1/M-expansion. It is simple to check that both types of additional term will lead 
to trivial or sub dominant contributions not relevant for our discussion. The Lagrangian 
displayed in (2) and (3) can be justified with an R-Symmetry, where both hidden sector 
fields X and Y have R charge |, iV has R charge |, E (rhd charged lepton superfield) 
has R charge 2 and the remaining superfields have R charge equal to 0, and the usual 
i?-parity is assumed, with in addition R P (X) = R P (Y) = +1. All dimensionless couplings 
g, h, etc, are taken to be order one parameters. 

Let us now suppose that after supersymmetry is broken in the hidden sector at scale 
raj the field Y acquires the following F and A-component vacuum expectation values, 

(Y) F = F Y = f Y m] 

(Y) A = A Y = 0, (4) 
while the field X acquires the F and A-component expectation values: 

\Xij)F = Fxij = 

(Xij) A = A X ij = axijmi. (5) 

Here f Y and axij are order one parameters, and the zero entries for Ay and Fx can be 
replaced by non-zero values Ay mj, and Fx <^ vrv\ without significant change. We 
will often re-write the scale toj in terms of the gravitino mass m 3 / 2 = mj/M, and the 
reduced Planck mass M. 

Before we proceed to analyze the consequences of the above effective Lagrangian for 
the light neutrino masses and mixings, a comment is in order concerning the assumption 
that (F x )ij mj. As is well-known, typical supergravity mediation of susy breaking 
has difficulties with FCNC and CP-violation constraints unless there is a high degree 
of degeneracy among the squark and slepton soft masses of different generations (we 
here ignore the possibility of alignment mechanisms which are typically much harder to 
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implement). Such degeneracy is accommodated in our model when (Fx)^ <C rnf. In this 
paper we will not be concerned about the detailed origin of this high-level of degeneracy, 
but take it as a phenomenologically necessary assumption. As we will argue below, in this 
case a simple prediction for the structure of the light neutrino mass matrix can result. 



3 Neutrino and sneutrino masses at tree level 

Given Eqs.(2) and (3), the relevant terms in the Lagrangian after the X and Y fields 
gain their expectation values are thus 

C = j d 2 6 (XijLiNjHu + M N NiN?) + AL^K + B^fi^ + ..., (6) 

where hi are the rhd sneutrino fields, Lj is the lhd slepton doublet, h u is the up-type 
Higgs scalar doublet, and the omitted terms include the usual soft scalar mass terms. 
First note that the effective neutrino Yukawa coupling in the superpotential of Eq. (6) is 
suppressed in magnitude by a factor of (m 3 / 2 /M) 1 / 2 ~ 10~ 7 — 1CT 8 , 



K=9axi j ]J 1 ^, (7) 

and gains its flavour structure from the (X^) A expectation value, and in addition, the 
scale of the rhd neutrino masses is lowered to the TeV scale 

M N = hf Y m 3/2 . (8) 

Second, there exists a TeV-scale, but flavour diagonal, trilinear scalar A-term 

A = g'f Y m 3/2 . (9) 

Finally there is a small but significant rhd sneutrino lepton-number violating B-term 
with coefficient 

Bl = h B f Y a Xij \[^ (10) 

of magnitude Bfj ~ (few x 100 MeV) 2 and flavour structure related to that of the neutrino 
Yukawa coupling {a 8 Xi ^ denotes the symmetric part of axij)- 

After electroweak symmetry breaking, with Higgs expectation values (H®) — v sin (5 
and (H®) = vcosfl, the effective Lagrangian Eq.(6) implies that at tree level in the 
neutrino sector, after integrating out the TeV-scale rhd neutrinos, there is generated a 
light neutrino mass matrix of see-saw type: 

M N 2 M 

If one takes M = 2 x 10 18 GeV, v = 174 GeV, and all a X ik, fy, 9 and ^'s to be of order 
1, this gives a contribution to the light neutrino mass matrix of order 10 _5 eV, certainly 
too small to generate the required Am^ and Amf 2 . 
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Note, however, that since the operators that give rise to the Yukawa and rhd Neu- 
trino mass terms are higher-dimension non-renormalizable operators, a more appropriate 
estimate of the couplings g and h might be the values found by applying the so-called 
'naive-dimensional analysis' (NDA) methodology which assumes that the cutoff M is 
bounded by the UV strong coupling scale of the non-renormalizable theory, and in which 
geometrical factors of l/(47r) 2 which enter loop calculations are taken into account[14]. 
(It is known that NDA works well for estimating the coefficients of the higher- dimensional 
operators in the low-energy chiral-Lagrangian description of QCD.) A simple calculation 
shows that the NDA estimates for the sizes of the couplings h and g which set the rhd 
neutrino mass and lhd-rhd-Higgs Yukawa coupling are h ~ 1 and g ~ 47r (and a UV 
strong coupling scale A related to the reduced Planck mass as A ~ 4irM) . This leads to 
an improved estimate of the size of the tree (see-saw) contribution to the light neutrino 
masses 

v 2 sin 2 3 g 2 

ii] — 



[m 



trcc\ 



10 _J eV. 



(12) 



M h 

Alternatively, the large mass scale M might be the string or susy-GUT scale, thus simi- 
larly increasing the estimate for m* ree . In any case, we will argue in the next section that 
a small tree-level term of size Eq.(12) can sometimes be an interesting perturbation to 
the dominant one-loop contribution to the light neutrino mass matrix. 

Turning to the sneutrino sector of the theory the effective Lagrangian Eq. (6) implies, 
after electroweak symmetry breaking, a 12 by 12 mass matrix that mixes the lhd and 
rhd sneutrinos and their conjugates (here, for simplicity, assuming A and B 2 real), 
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where we work in the basis (v*,n, v,n*), Ml = m\ + m|cos(2/3)/2 is the lhd sneutrino 
mass arising from the usual soft mass m\ and electroweak breaking terms, and M R = 
m 2 R + Mlj is the total rhd sneutrino mass including the soft mass-squared. (Our notation 
is that Latin indices run from 1 to 3, so e.g., the n are indexed as % + 3, and sneutrino 
mass eigenstates are labeled by Greek indices a, (3 = 1, ... ,12. ) This mass matrix is 
diagonalised by a unitary rotation of the form 



U = 



( 



V 



V/y/2 


-V/y/2 






V/V2 


-V/y/2 



V/y/2 


V/y/2 




\ 

V/y/2 


V/y/2 ) 



( 



V 



COS</) 

- sin i 





sin0 

cos(/> 

COS0 







sin 

sin0 cos0 J 



(14) 



where tan 20 = 2Av sin /3/(M| — M R ), and V is the 3 by 3 matrix that diagonalises 
B 2 j. The expression Eq.(14) is correct up to small terms of order B 2 j(M\ — M R ). The 
corresponding sneutrino mass eigenvalues fall into 2 groups of almost 3-fold degenerate 
complex states. 



5 



Bip 

-x - - - - 



A8 kk> / 



•xA8 m 



V 



Figure 1: The dominant contribution to m u . 



4 Structure of light neutrino masses at 1-loop 

The rhd sneutrino interaction Eq.(lO) gives rise to a radiative contribution to the light 
neutrino masses, illustrated in Fig. 1, that, for (Y) F ~ m 2 and (Xij) A ~ mi, dominates 
over the tree level contribution arising from the the mixing of the TeV-mass rhd neutrinos 
with the lhd neutrino states. In detail, the contribution to the light Majorana neutrino 
mass matrix from the diagram in Fig. 1 is given by, 

= ~ iq 2 ~^a,*+6^1,jH-6^/3,AH-3^/!^ P, 1, 8, x), 

(15) 

where all repeated indices are summed over. In this expression, Xx is a factor that depends 
on the exchanged neutralino, U is the sneutrino mixing matrix, x — 1, ...,4 denote the 
neutralino mass eigenstates, and L(a, j3, 7, 5, x), is a totally symmetric function of the 
sneutrino and neutralino masses that arises from the momentum integral in Figure 1. 

Because of the small effective Yukawa coupling, Eq.(7), between the higgsinos, lhd 
neutrinos, and rhd sneutrinos, in practice only the bino and zino components of the 
neutralino are exchanged across the bottom of the loop. Thus only these states need to 
be expanded in terms of the 4 neutralino mass eigenstates, leading to 



2Mf 



v 2 



N*i sin 9 W - N* 2 cos 9 W 



M Xx (16) 



where N xI is the standard neutralino mixing matrix, and M Xx are the neutralino masses. 

A transparent and elegant form for the light neutrino mass matrix results if we 
consider the situation in which the A-term mixing is small. This greatly simplifies the 
flavor structure of the result Eq.(15). (The V matrices in U of Eq.(14) cancel by virtue 
of VV^ = 1.) Moreover, in this limit the loop factor L also simplifies, 

1 = Wr - mr- mi - *Y [ (1 " r>(1 + r - 2x){1 ~ x){ - ~ x) + ( 2f2 " xHl + r> 

+ 2x 2 (l + r 2 + r^J logr - a;^log(r/a;) +rlog(x 3 r) - r 2 \og(x 3 / r A ) +r 3 loga;^ , (17) 
where r = M\jM\ and x = M^jM\. Therefore in this limit of small sneutrino mixing 
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the 1-loop contribution to the light neutrino mass matrix becomes 



TO, 



Xx A 2 B 2 jV 2 sm 2 p 



)« = E - — — M i M ^ M l> ( 18 ) 



The most important feature of the result Eq.(18) is that the overall scale of the 
contribution m l ° op is naturally of the correct size to account for atmospheric neutrino 
oscillations (as is also true of the models of Refs.[8, 10]). To see this explicitly it is useful 
to consider the simple case in which all the lhd and rhd sneutrinos and the neutralinos 
are approximately equal to a common mass scale, m susy , giving 

(19) 



12< usy 



which leads to a one-loop contribution of magnitude 

m'r = - 10 ~ 2 eV - 10-1 eV (20) 

967T M b m° nsy 

depending on the precise magnitude of the A and B terms. In addition, in our model, 
the flavour structure of this dominant one-loop contribution to the light neutrino mass 
matrix is determined directly and entirely by the rhd sneutrino lepton-number violating 
S-term, B 2 j, which is in turn generated by the {Xij) A expectation value. Moreover, as 
claimed earlier, the one loop result dominates over the tree level see-saw contribution. 
It is useful to define the (small) parameter e as the ratio of magnitudes of the tree-level 
see-saw contribution Eq.(12) to the above 1-loop contribution. 

It is also interesting to consider the regime in which the rhd sneutrino states are 
heavy compared to the neutralino and lhd sneutrino states, r ~ x <C 1. In this case the 
loop factor is approximated by the expression 

L (x-r)+xlog(r/x) 1 

M R (r-x) 2 ' 1 ' 

Taking Ml ~ M x ~ m susy and scaling the A and B terms relative to their natural values, 
Al ~ m 2 usy , and B% ~ m 2 usy (m susy /M) 1/2 given by Eqs.(9) and (10), leads to 



167r 2 M|MV2 \A 2 

Assuming m susy ~ 300 GeV, and M R ~ 1 TeV this gives m l ° op ~ 0.02 eV(A 2 /A 2 )(B 2 / B%) 
showing that M R cannot be much heavier than 1 TeV unless the scale of the MSSM 
superpartners is uncomfortably high. 

In either case the final structure of the light neutrino mass matrix is in total 

K°%- = fi(a s x + ea T x a x )ij. (23) 

with the scale set by \x ~ 0.1 eV — 0.01 eV, and e in the range e ~ 10~ 2 — 10~ 4 . An 
attractive feature of this structure is that it allows us in a simple way to account for 
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the hierarchy between the solar and atmospheric neutrino mass-squared splittings. The 
atmospheric Am 2 can arise from the one-loop contribution, while the tree-level correc- 
tion leads to the small Am 2 olar splitting, the hierarchy being entirely due to the small 
dynamical parameter e ~ 1CT 2 . Of course it is possible that the hierarchy instead arises 
entirely from the flavour structure of the dominant [xa x term, the tree-level perturbation 
being insignificant, but this leads to quite traditional neutrino flavour models, so we here 
focus upon the new possibility. 

As a simple example of a model along these lines consider the situation in which the 
matrix a x corresponds to the leading-order form for an inverted hierarchy model 



a x 




(24) 



Therefore in this case the full light neutrino mass matrix including both loop and tree 
contributions has the form 




(25) 



This mass matrix has one zero eigenvalue, and two massive eigenvalues m± 
n(\/a 2 + 1 ± e). Therefore to accommodate the oscillation data requires (1 + a 2 )// 2 



2 x 10 3 eV 2 , while the solar oscillation data requires 4(1 
a ~ 1 this gives 



■a 2 ) 1 / 2 /^- 



/i 2 ~l(T 3 eV 2 , e~l(T 2 



7x 10~ b eV z . For 



(26) 



comfortably of the sizes expected in this model. Moreover, since m* ot is diagonalised by 
first performing a 23 rotation, i?23(#), with angle 9 = tan -1 (a), and then a 12 rotation, 
Ri2(<f)), with = it/ 4, we see that for a ~ 1 the atmospheric and solar angles will both be 
close to maximal. Specifically, recalling that the physical neutrino mixing (MNS) matrix 
includes a unitary matrix Vl from the rotation of the charged leptons to a basis where 
the charged lepton mass matrix is diagonal and real we find a MNS matrix of the form 



^mns = V L Rl z {9)RU<t>) = Vl 



( S0 \ 

-cqs^ cqc^ s e ■ (27) 

\ SgS^ -SgC^ Cg J 

Under the not unreasonable assumption that the mixing angles in from the charged 
lepton sector are small, Vl only slightly perturbs the above structure, thus leading to 
almost maximal physical atmospheric and solar mixing angles, and a small 6> e3 angle. 

Alternatively, light neutrino masses with a normal hierarchy can be easily generated 
if a x contains an antisymmetric piece, a x . In particular consider the forms, 
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Substitution of these matrices into Eq.(23) with reasonable (and not fine-tuned) 0(1) 
values for the parameters b, c, d, and e, and values of the dynamically determined pa- 
rameters /i 2 ~ few x 10" 3 eV 2 , and e ~ 1CT 2 , leads to a light neutrino mass matrix, which 
when diagonalised produces mass squared differences within the experimental bounds for 
a normal hierarchy. Moreover, again under the assumption that the real mixing angles 
in Vl are small this leads to large physical atmospheric and solar mixing angles, and a 
small # e3 angle. More generally, if the matrix ax is not real as in the above examples, 
but contains large phases then it is simple to generate successful models in which the 
one-loop term gives Am 2 tm while the perturbing tree term leads to Am 2 olar . 

5 Comments and conclusions 

In this paper we have further analyzed a class of models first introduced by Arkani-Hamed 
et al [7, 8], and Borzumati et al [9, 10], in which the light neutrino masses are a result of 
higher-dimensional supersymmetry-breaking terms. The mechanism is closely related to 
the Giudice-Masiero mechanism for the MSSM Higgs /i parameter, and in particular leads 
to TeV-scale rhd neutrino and sneutrino states, that are in principle accessible to direct 
experimental study, unlike traditional see-saw mechanisms. A second difference is that 
the dominant contribution to the light neutrino mass matrix (which is of Majorana type) 
is a one-loop term induced by a lepton-number violating and supersymmetry breaking 
B-term for the sneutrino states that is naturally present in the model. In this letter 
we have focused upon the simplification and analysis of the flavour structure of this 
general class of models, and have found that simple predictions for the light neutrino 
mass matrix are possible. In addition we have found that the subdominant tree-level 'see- 
saw' contribution may lead to interesting perturbations of the leading one-loop-induced 
flavour structure, possibly generating the smaller Am 2 olar . 

In this paper we have not explored the important issues of the possible collider and 
cosmological tests of our models. In broad structure the implications of our models are 
similar to those already analyzed in Refs.[7, 8, 9, 10]. In particular one expects the A- 
term interactions in our model to lead to interesting possibilities for production and 
decay of the TeV-scale rhd sneutrino states, including the possibility of anomalous Higgs 
decays. Another intriguing possibility is that the rhd sneutrino states could be the dark 
matter [15], or that, when CP- violation in the neutrino sector is taken into account, 
dark matter with dominant inelastic interactions with matter results [16]. In a future 
publication, [17], we show that the class of models analyzed here naturally lead to a very 
attractive and successful theory of TeV-scale resonant leptogenesis, developing from the 
earlier work of [18]. 
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